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Q ■ 1. Introduction 

One of the main tasks of the quantum statistic physics is calculation of the correlation functions 
for many-body systems of different kind because they contain the most important information 
about the observable quantities and system properties. As it was first noticed by Kubo [Q], linear 
^SJ ' transport coefficients are expressed in terms of the Fourier transforms of appropriate correlation 

^O I functions, which relate by spectral relations to the two-time Green functions. Since that time 

^T I the Green's function method has been noticed and extensively developed [EIEIS], that provides 

equations which allow to calculate the expectation values of operators and observable quantities. 
("^ I But very soon it was noticed that spectral relations must be completed by the special treatment 

of the pole at zero frequency, which gives additional contribution connected with presence of the 
C^ ' conserving quantities [EIIEIIZI- Later, it was shown that such contributions describe the difference 

between the isothermal and isolated response of the many-body system and are specific for the 
non-ergodic systems [^IBJOIS] where the regions in the phase space which can not be achieved by 
the trajectory of the point that describes an evolution of many-body system exist. In the Green's 
function formalism the issue of ergodicity appears as a difficulty in the determination of the zero- 
O '■ frequency bosonic propagators [|niini|7||ni[TOl[TTlIiaiISIIIlIiniIiniiniIIH|- Nevertheless, even 

now many textbooks on the quantum statistics and many-body theory do not provide complete 
discussion of the spectral relations and special treatment of the zero-frequency functions. 

The equations of motion do not uniquely determine the causal and retarded Green's functions. 



X 

^ . but only up to J-function of frequency with some unknown coefficients which produce additional 

contributions in the zero- frequency functions (see, section |21). Usually, these zero-frequency func- 
tions are fixed by assigning them their ergodic values, but this cannot be justified a priori. A 
wrong determination of them dramatically affects the values of directly measurable quantities like 
compressibility, specific heat, and magnetic susceptibility. In order to handle this zero-frequency 
functions a different approaches were developed, e.g. anti-commutator bosonic Green functions [ 
I14| . direct algebraic method [E|) singular- value decomposition [^3, algebra constraints in the 
composite operator method [121122) etc. 

On the other hand, temperature Green functions [ I23L ^^ are free from this issues and allow 
to avoid all complications connected with the presense of the non-ergodic terms [Eji e.g. they 
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allow to calculate isothermal susceptibilities for the Ising model, where Kubo response is equal to 
zero, as well as for the more complicated spin and electron [ I26LI27J and pseudospin-electron [ I28| 
systems, and for the infinite-dimensional Falicov-Kimball model [|551|2n|j where exact expression 
for the isothermal charge susceptibility [' 31[I32| contains both Kubo response, which is finite at all 
temperatures, and non-ergodic contribution, whose divergencies give the phase transition points. 
The mentioned above issues about special treatment of the zero-frequency functions are mostly 
investigated for the two-time correlation and Green functions, but nobody have cosidered this for 
the multitime one. In the original Kubo's paper [HI, the formulation of the transport theory is not 
limited to the linear phenomena and solution of the Liouville equation for the density matrix is given 
to the arbitrary order in the strength of disturbance. Resulting multitime correlation and Green 
functions can be used for the description of the nonlinear transport phenomena and resonances [ 
I33L 1341 IS5] . Besides, multitime correlation functions also appear as puzzles in different orders of 
the perturbation theories for many-body systems [131]. Moreover, cross-sections of the inelastic 
scattering processes can be expressed in terms of the multitime correlation functions too, e.g. for 
the electronic inelastic light (Raman) scattering the nonrcsonant, mixed, and resonant responses 
are connected with the two-time, three-time, and four-time temperature Green functions [ I37[I38| . 
respectively, and can be rewritten in terms of the multitime correlation functions. Nonresonant 
contribution is connected with the spectral density of the two-time correlation function 

RN{q, n) - 2Trg\k,)g\ko)I^4n, ~n), (1.1) 

mixed contribution is connected with the spectral density of the three-time correlation functions 



RMiqM) = 2ng^ih)g^iko) f duj 



+ 00 

'^^,(o),(i)(f^,CJ-^, -W) 



(1.2) 



uj + ujo — iS to — uji — iS uj + LUo + '^S 

and resonant contribution is connected with the spectral density of the four-time one 

RR{q,n) = 2V(fc,)g2(fc„) f diU f dQ (1.3) 

— CXD — OO 

r Jj(i)j(o)^(o)j-(,) {u,fl- u;,u} - ri, -(!)) Ij(o)J(^)J(^)J(o) {uj , fl - UJ , u} - fl, -u}) 

{lO ~ LUi — i(5)(cZ) — UJi + 16) (lU + LUo + iS){LL! + LUo — i(5) 

Ij('')j(i)j(o)j(^) iuj,fl- uj,u} - fl, -u}) /j(i)j(o)j(>)j(o) (cj, il - tj, a) - ri, -(!)) 



{u + L1J0 + iS){u; — UJi + id) (uj — uJi — iS){uj + uJo — iS) 

where /yiB...(u;i,a;2, . . .) are spectral densities of multitime correlation functions (see below) and 7 
and j^^'°' are expressed in terms of the stress tensor and current operator, respectively. 

In general, spectral relations include (a) representation of the observable quantities in terms of 
the spectral densities of the multitime correlation functions [like Eqs. (|l.l(l . (|1.2|l . and H1.3(l for the 
inelastic light scattering], (b) extraction of the non-ergodic contributions in the spectral densities, 
(c) representation of the multitime temperature Green functions in terms of the spectral densities, 
and (d) solution of the reverse problem — finding of the spectral densities from the known multitime 
temperature Green functions and calculation on this basis of the observable quantities. Spectral 
relations for multitime correlation functions were formulated in Ref. [ 01 |3HI only for the case of 
the multitime generalizations of the retarded and advanced functions without any consideration 
of the zero- frequency non-ergodic contributions. The main purpose of this paper is to show how 
non-ergodic contributions enter in multitime correlation functions and how a complete set of the 
spectral relations for the multitime correlation and Green functions can be derived. 
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2. Two-time correlation functions 

Before we start consideration of the spectral relations for the multitime correlation functions, 
let us remind basic relations and main results for the two-time correlation functions. In general, 
two-time correlation function is defined by the following expression 

KABih~t2)^{A{h)B{t2)), (2.1) 

where A{t) = c'^*Ae~'^* is an operator A in the Heisenberg representation, and angular brackets 
denote statistical average with total Hamiltonian of the many-body system. Here and below we 
shall consider the case of the bosonic operators only, generalization for the fermionic one is obvious. 
We can perform Fourier transformation of Eq. (|2.1|) that gives a spectral density for the two-time 
correlation function 

Iab{uji, -uji) = ^ fd{h - t2)e''^''^''-''^KAB{ti - ^2) (2.2) 

— 00 

= lAsii^l, ~^l) + (5(wi)/ab(o, o), 

where we have separated two contributions 

iABiu^i,-iUi) = ^^'e-''^^A,/B,,,5(e,, -1-c^i), (2.3) 

Iab{o,o) ^\Y. e^'^'^'^J'^y (2-4) 

with different frequency and time dependences. The first one is time-dependent and includes sum 
over states with different energies [eji = Ej — ei ^0) denoted by prime 1)2. 3|l . we shall call it 
regular contribution, and the second one is purely static and includes sum over the states with 
the same energy H2.4f) . we shall call it non-ergodic one. Here, Aji = {j\A\l) are matrix elements 
of the operator A, Z = ^j e~^'^' is partition function, and o denotes that given contribution 
does not depend on the respective frequency. Besides, spectral densities (j2.2|l satisfy the following 
permutation relation 

lAB{i^u-toi) = lBA{-tui,LOi)eP^' (2.5) 

and for the given operators A and B there is only one nonidentical spectral density. At this point 
there are no special need to separate regular and non-ergodic contributions, and, as a rule, in 
the textbooks nobody do this. But such separation becomes important when one is going to find 
spectral densities using spectral relations for the two-time retarded Green's function 

G^t(ii - h) = {{A{h)\B{h))) = -■Mh-h){[A{t,),B{t2)]). (2.6) 

Its Fourier transformation, from the formal point of view, 

oc 

Gab{coi, ~LOi) = ^ /"d(ti - t2)e'^'^''-'-'^GAB{ti - t2) (2.7) 

= Gab{(^i, -ioi) + (5(wi)Gc(o, o) 

also includes two contributions with different time dependences 

+00 
GABiiOi,-uji) ^ / da)i/As(tZ'i,-a)i)-: -— , (2.8) 

J LJl — (jJl ± 10 

— C30 

5^5(0,0) =/^b(o,o), (2.9) 
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but only first one can be derived using an equation of motion technics 

ojiGAsii^i, -^i) = ^iGab(wi, -uji) (2-10) 

and the second one is omitted in the equations of motion 

U!iS{uJi)GABi°,o) = 0. 

As a result, a corresponding static contribution into spectral density can not be handled directly 
by the spectral theorem for the retarded Green functions. 

On the other hand, a method of the temperature or Matsubara Green functions allows to avoid 
all complications connected with the presense of the non-ergodic terms and obtain all contributions 
in the spectral density in a straightforward way [jlSj. In general, a two-time temperature Green's 
function can be defined as 

Kcin - T2) - (ri(ri)B(r2)), (2.11) 

where A{t) = e^'^Ae~^'^, r is imaginary time (inverse temperature), and T is operator of the 
imaginary time chronological ordering. Its Fourier transform 

P 
K,(ii/i, -iiyi) = y"d(Ti - T2)e''''(^'-^''>K,{n - r^) (2.12) 


= kcCwi, -ivi) + /3A(iz/i)^e(o, o), 

where \vi = ilo,^-^ = 2t:[Ti>i are Matsubara's frequency and 

A(.)^{J; 1 = 1 (2.13) 

is generalization of the Kronecker symbol, also contains two contributions with different time 
dependences 

Kc{ivi,-\vi) = / du)ilAB{uJi,-u)i)— : , (2.14) 

J LUl- lUl 

— C30 

i^c(o,o) = /^s(o,o). (2.15) 

As a rule, temperature Green functions are handled using different kinds of the diagrammatic 
technics which allow to calculate contributions with A-symbols and according to (|2.15(l they can 
be identified as a non-ergodic contributions in spectral densities, that is the main difference with 
the case of the retarded Green functions. For the non-ergodic systems such contributions are the 
main one which determine critical behavior. A typical example is the Ising model for which there 
are only non-ergodic contributions whose divergencies give Curie points and there are no regular 
one. Another example of the non-ergodic fermionic system is the Falicov-Kimball model [ 1291 130| 
whose isothermal charge susceptibility contains both non-ergodic and regular contributions but 
only the first one determines the critical point [ 1311 132 | . Regular contribution in spectral density 
is connected with nonanalyticy (imaginary part) of the Green's function at real axis and can 
be obtained by performing an analytic continuation of the Matsubara Green's function from the 
imaginary to complex frequencies and then to real one iz^i — > zi — > wi ± i(5 



= ^^=(^i'~^i) 



^1— C(Ji+i(5 



= /AB(u;i,-c^i)(l-e-'3-i), (2.16) 

2i— cji— i(5 



Kc{zi,-zi) 
that completes a spectral relations for the two-time correlation functions 
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3. Three-time correlation functions 

Now let us proceed to the consideration of the three-time correlation functions. Generalization 
for the case of the higher-order multitime correlation functions can be done in the same way, but 
it is much more cumbersome and will not be considered here. 

Three-time correlation function can be defined in a usual way as 

KABcitiMM) = {A{h)B{t2)C{h)). (3.1) 

Here, we shall consider only the case of the equilibrium many-body systems which are time-shift 
invariant 

KABc{ti,h, h) = KABcih -tM- t, h ~ t). (3.2) 

Spectral density is defined as its Fourier transform 



(^i,c.2,c^3) = 7^ j d(ti -is) j d{h-h)e'^'^^*'+'^-'-+'^-'^^^KABc{tiMM) (3.3) 



= [/aBc(wi,W2,W3) +(5(wi)/asc(o,-W3,W3) + 5{Uil2)lABc{^l,°-, -^l) 
+ 5(cJ3)/abc(-W2,CJ2,o) + 5(cJi)(5(w2)iASc(°,°,o)]A(a;i +UJ2 +0^3) 

and includes five different contributions with different time dependences 

/abc(wi,w2,-wi -W2) = — ^ e-^^^ AjiBifCfj5{eji + uji)5{eif + UJ2). (3.4) 



Iabc{o.-ojz,^z) = - Y^ e l^'^AjiBifCfjSiefj+ujs), (3.5) 

jif 

iABc{toi,o,-uj^) = ^ Y^ e-'^'^AjiBifCfjSieji+coi), (3.6) 

Iabc{-^2,oj2,o) = - Y (^~'^'''^jiBifCfjS{eif+uj2), (3.7) 

^ jif 

/ABc(°,o,o) = i Y ^'^"'AjiBifCf,. (3.8) 

Besides, the total spectral density (|3.3|l as well as each contribution satisfy the following cyclic 
permutation identities [uji -I- 072 + ^3 = 0) [flll36| 

lABc{<^l,^2,^i) = /BCA(w2,W3,Wi)e'^"i = /cAb(w3, ^1, ^2)6"^"^ (3.9) 

and for the given operators A, B, and C there are only two nonidentical spectral densities, e.g. 

lABc{^l,^2,i^3) and IcBA{i^3,i^2,i^l)- 

Now we introduce three-time temperature Green's function 

^c(Tl,r2,T3) = (ri(ri)B(T2)(7(T3)), (3.f0) 

Kc{ti,T2,Tz) = Kc{ti - T, T2 - T, Ti - t) . 

Due to the imaginary time ordering its Fourier transform contains 3! = 6 terms which can be 
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collected into two groups of three terms connected by the cyclic permutations 

f3 13 

K,(i^i, ii/2, i^^a) = ^ / d^i / d^2 1 dT3e('''i^^+"^=^^+'''^^^)i^,(Ti, T2, ra) 



(3.11) 







2 X! [^i'^'/'^/J*P(-?' i^i' ^' i^2, /, ii^s) + CjfBfiAij^iJ, iz/g, /, iL'2, 1, it^i)] , 



jU 



where 



P Tl T2 

— R^- I , 



(3.12) 



/3 T3 T3 



^"^ /dn /drs f dT3+e-^'' f dT2 /"drg /" dn + e"^^^ /"drg /"dn /dra 





X exp[(ej7 + iiyi)Ti + {eif + 11^2)^2 + (e/j + ii^3)T3]- 







In the general case, when all Matsubara frequencies are nonzero or when there are no eigenstates 
with the same energy value, function H3.12|l is equal 



= A(ii^i + W2 + ii^s) 



(3.13) 



-/3ej 



-Pei 



-pes 



+ 



+ 



{eij- 'ivi ) (e/j + ij^3 ) (e/i - i J^2 ) {sji + wi ) (e»/ - W3 ) {eif + ii/2 ) 



Besides, we must consider several special cases, when we have levels with the same energy value: 
case £j = ei ^ Ef and wi — —\V2 — ii^s = 0, when 



<P(j, 0, 1, W2. /, -yv2) - A(ii^2 + lyz) 
case El — Ef ^ Ej and ii^2 — ^11^3 — iiyi = 0, when 

*P(i, -i'^3,^,0,/,iz^3) = A(ii/3 + ij/i) 
case £/ = Ej 7^ £; and ii^a = — ii/i — ii^2 = 0, when 

*P(j,ii^i,^,-ii^i,/, 0) = A(ii^i +ii/2) 
and case e-, = e; = £/ and ii^i = ii^2 = ii^3 = 0, when 



(3e- 



^Siei 



'H^f 



-Pei 



£fi ~ ii^2 (e/i - iJ^2)^ 

— \ 



Elj - \Vi [eij - WiY 



<P(j-0,Z,0,/,0) = ^e-'5^^. 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



The second term in the r.h.s. of Eq. H3.14|l can be derived from Eq. H3.13|l by the analytic continua- 
tion of the Matsubara frequencies to the complex one iz^i -^ zi and it/3 -^ — ii^2 ~ Zi followed by the 
limit zi ^ 0, but the first one can not be derived from Eq. (|3.13|) by manipulating by frequencies 
only and corresponds to the additional non-ergodic or conserving contribution, which originates 
from the presence of the states with the same energies, and appears only at zero frequency. Such 
non-ergodic contributions can be derived from Eq. H3.13|l by 



lim lim — lim lim 

Eji— >0 2i^0 2i— ^OEj;— *0 



^{j,Zl,l,Z2,f,-Zl - Z2) = 



(3e- 



-I3ei 



Efl - Z2 ' 



(3.18) 



but such derivation involves manipulations with the many-body quantum states energies, that can 
not be, in general, reproduced by the quantum statistics many-body methods. A similar analysis 
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can be done for Eqs. H3.15|l and H3.16|l and, after analytic continuation from the imaginary axis to 
the complex plane ivi —^ Zi with using constraint 



^^^ = 0, 



(3.19) 



Eq. H3.12(l can be rewritten as 



^(i, Zi,l, Z2, f, Z3) = A(zi + Z2+ Z3) 



^A(zi)A(z2)A,^,,„,,e-''^^ 



-Pcf 



(3.20) 



-l3ei 



+ /3A(zi)Ae^.,,, ^ + /3A(z2)Ae,,e, 



^ Sfl - Z2 



.-Pej 



-I3ei 



Q-Pej 

— + /5A(z3)Ae^,,^ 

eof-z^ 'eij-zi 

Q-PSf 



{eij- zi){efj+z:i) {efi-z2){e,i + zi) {e^f - Z3){eif + Z2) 



where 



^Ej, ...,£/ 



1, ej=...=ef 
0, other case 



(3.21) 



Now, after substitution of Eq. (|3.20() in Eq. (|3.11l) . we get the following representation for the 
analytically continued three-time temperature Green's function 

K,{zi,Z2, Z3) = ^A(zi)A(z2)A(^3)i^c(o, o, o) + /3A(zi)A(^2 + ^3)1^0(0, Z2, Z3) (3.22) 

+ l3A{z2)A{z3 + zi)Kcizi,o, Z3) + /3A(z3)A(zi + 22)^0(^1, ^2, o) 

+ A{zi + Z2+ Z3)Kc{zi,Z2, Z3). 

It includes five contributions which always can be distinguished by the different A factors. The 
first contribution is expressed directly by the spectral densities of the H3.8|l type 

Kc{°, o, o) = iABc{°, o, o) + icBAi°, °, °), (3.23) 

the next three contributions are expressed in terms of the spectral densities of the (|3.5|I - H3.8|) type 



^c(o,^2,-22) 



da;2 



iABc{°, X2,-X2) - icBA{~X2,X2,o) 



X2 — Z2 



(3.24) 



— [-^Asc (0,0,0) --^"05^(0,0,0)], 

Z2 



+00 



Kc{-Z3,o,Z3) 



da;.' 



Iabc{-X3,°, a;3)eP^3 _ Icba{x3, o, -X3)e ^^^ 



PX3 



X3 - Z3 



(3.25) 



and 



H [iABc{°, o, o) - icBA{°, o, o)] , 

^3 



+00 



Kc{zi,-zi,o) = / dxi 

— C 

1 



i ABC {:Xl,-Xi,o) - icBA{°, -Xi.Xi) 



Z\ 



X\ — Z\ 
[iASc(o,o,o) -/cba(o,o,o)], 



(3.26) 
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and the last contribution is expressed in terms of the spectral densities of the H^14|l - (|^17|l type 

Iabc{o, X2, -a;2)(l - e-^^-) + /csa(-X2, 2:2, o)(l - e''-^) 



Kc{zi,Z2,Z'i) = / da;2 



+ 00 



dX3 



lABC 



(Z2 - X2){zz +X2) 

(-a;3,o,x3)(e''-^ - 1) + /cBA(a;3, o, -X3)(c-''--^ - 1) 



(3.27) 



—00 
+00 



dxi 



(Z3 -X'i){zi +X'i) 

iABc{.xi,-xi,o){l - e-l^-^) + icBA{o,-xi,xi){\ - el^^^) 



— C50 

+00 +00 



dx3 / dxi 



— 00 —00 

+ 00 +CXD 



dxi / Ax 



—00 —00 
+00 +00 



da;2 / dx 



(Zi -Xi){z2 +Xi) 

lABcixi,-X3 -Xi,X3) + IcBAix3,-X3 -a;i,Xi) 
(2^3 " Z3){xi -Zi) 

iABcixi,X2,-xi -a;2)e~^^i + Icba{-xi - X2,X2,xi)e'^''^ 

(xi - Zi){x2 - Z2) 

iABc{-X2 - X3,X2,X3)e'^''-' + JcbaJxs, X2, -xi - a:2)e~^^3 

{X2 - Z2)ix3 - Z3) 



Here, we have used the cyclic permutation identities (|3.9|l according to which there are only two non- 
identical three-time correlation functions, e.g. Iabc{^i^'^2,i^3) and IcBAi^3,^2,^i)- Eqs. (|3.22|) - 
H3.27|) give complete representation of the three-time temperature Green functions in terms of the 
spectral densities. 

Now we pass to the solution of the reverse problem — finding of the spectral densities from 
the known multitime temperature Green functions. According to Eqs. H3.22|l - H3.27|l each of five 
contributions in the three-time Green's function (|3.22|l has different set of branch cuts, that allows 
to extract all ten spectral densities, which enter. 

First of all we perform an analytic continuation of Eq. H3.24|l (z2 — > ^2 ± 16, Z3 — —22, zi =0): 



Kc{o,Z2,-Z2) 



= -Kc{o,-Z3,Z3) 



= S{UJ2)S^{0,0,0) +Jii{o,UJ2,-LL!2), 



(3.28) 



where 



i^(o, o, o) = lABcio, o, o) - icBA{°, o, o), (3.29) 

^l{o^UJ2,-UJ2) = /aBc(o, ^2,-^2) - icBAi-UJ2,t^2,°)- (3.30) 

Next, we perform an analytic continuation of Eq. I|3.25|l (23 — > cl;3 ± iS, zi = —2:3, Z2 = 0) 



Kc{~Z3,o,Z3) 



= -Kc{zi,o,-zi) 



= -(5(cj3)i^o (0,0,0) +il2(-t^3,o,Ci;3), (3.31) 



J?2(-'^3,o,CJ3) = /ABc(-'^3,o,CJ3)e'^"3 -/CSA(W3,0,-W3)e ''"^ 

and of Eq. (|3.26|l (zi — » wi ± iS, Z2 = — zi, 23 = 0) 



(3.32) 



Kc{zi,-zi,o) 



= -Kc{-Z2,Z2,^ 



^ S{uji)Mq{o, o, o) + il3(a;i, -wi, o), (3.33) 

il3(a;i,-wi,o) = /^Bc('^i,-'^i,o) ~ icBA{°,-^i,(^i)- (3.34) 

One can see, that in Eqs. (|3.28() . H3.31|l . and H3.33|l all factors at (5-functions are the same. 
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Next step is more complicated and requires a two-stage procedure. First of all we perform an 
analytic continuation of Eq. (|3.27|) over first frequency zi —^ uji ± iS (za = — wi — Z2) 



Kc{zi,Z2,Z3) 



= ^ [lABci^i, o, -uJi){c-^^' - 1) + /csa(-^i, °, wi)(e''"^ - 1)] (3.35) 

■ [lABciuJu -wi, o){e-P^' - 1) + icBA{o, -L^i,c^i)(e^"^ - 1)] 

iABc{^l,X2.-UJl - .T2)(e-''"i - l)+icBA{-OJl - .T2, ^2, ^i) (e'^'^^ - 1) 



Z2 

Z2 +^1 

+ OC 



dx2- 



Z2 -X2 



and then an analytic continuation over second one Z2 ^ 0^2 ± i5 (23 — > —uji — 0^2 T i<^) 



Xc(zi, 22,^3) 



-Kc{zi,Z2,Z3) 



;(2) 



(3.36) 



(3), 



(5(cj2)i^iJ(-t^3,o-^3) + 5{u}3)^ia{^i,~^i,°) + i?i, 2(^1,^2,^3) 



where 



•'?o(-^3,o,W3) 



-/ABC(-W3,0,^3)(e'''"^ - 1) -/cBA(t^3,o,-W3)(e~'''"^ " 1), 

R^liuji, -c^i, o) ^ ^lABcii^i, ~u;i,o){e-^'^' - 1) - icBAio, -c^i,o.i)(e'5'^i - 1), 

J^l,2(wi,ti^2,'^3) = -/ABc(wi,'^2,W3)(e"^'^i - 1) - /cBa(w3, t^2, Wi)(e'^"i - 1). 



(3.37) 
(3.38) 
(3.39) 



We can perform another sequence of the analytic continuations: Z2 — > a;2 ± iS (zi = —072 — Z3) and 
Z3 -^ UJ3± 16 (zi -^ ~UJ2 — W3 =F i(5), that gives a different set of relations 



Kcizi,Z2,Z3) 



= -Kc{zi,Z2,Z3) 



= 5{U!3)^1^{UJI, ~LUl,o) +5{uJl)^l{o,UJ2,-L02) + J?2,3 (^1 , CJ2, W3) , 



(3.40) 



where 

4^^(c^i, -c^i, o) = -lABciuJi, -oJi, o)(l - e-'^'^O - icBA{o, ~uJi,0Ji)(l - e^"0, (3-41) 

.^^) (o, U2, ~L02) = -lABcio, ^2, -^2){e-^'^' - 1) - /csa(-c^2, ^^2, o)(e'3-^ - 1), (3.42) 

.^2,3(^1,^2,6^3) = -/Asc(c^i,^2,W3)(e'^"^ -e-''"0-^CBA(^3,^2,c^i)(e"'^"^ -e'3'^i). (3.43) 

One can see that 

^ Kc{zi,Z2,Z3) 



Kc{zi,Z2,Z3) 

but there are common elements in final expressions, e.g. 

^1^1(^11-^1,°) +-*?2^3(^l'^^l'°) = 0, 

and the following Jacobi type identity is fuUfilled 



ii:c(zi,z2,z3; 



Kcizi,Z2,Z3] 



+ Kcizi, Z2, Z3) 



0. 



(3.44) 



(3.45) 



(3.46) 



Due to this identity the procedure is unambiguous and resulting set of equations for the spectral 
densities is not overdetermined. 

In general, the three-time temperature Green's functions H3.22|l . as well as multitime Green 
functions of higher order, have very complicated analytic properties as function of the complex 
frequencies Zi, but they can be separated into contributions with different A- factors and different 
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frequency dependences. Analytic properties of this contributions can be also complicated, never- 
theless, we can always derive a complete set of equations for the multitime spectral densities from 
the multitime temperature Green functions by the different sequences of the analytic continuations 
Zi ^ LUi zt iS accompanied by the constraint (|3.19() . 
Finally, we find from Eqs. (pT^ and ^2^: 

lABcio, °, °) = 2 [^-(°' °' °) + M°, °, °)] , (3.47) 



Icba{o, °> °) = 2 t^'^^"' °' °) ~ ^^°' °' °)] ' 



from Eqs. ifX^ and (|^^ : 

Icba{-'UJ2,uj2,o)= ^p^^_^ + ■ (e/3... _i)2 ' 

from Eqs. l|0^ and (|07jl : 

f . ^ ^2(-W3,0,W3) , 4^^(-W3,0,W3) 

/ABC (-^3, 0,^3)= ^^^^_^ + ^^^^3_^^, , (3.49) 

? , Jl2(-^3,°,W3) , ■^1,2(-'^3,°,W3) 

1cba{^3: °, -'-'3) - 1 _ e-/3^3 + (i_e-/3-3)2 ' 
from Eqs. ifX^ and (|^3S|) or jX^ : 

JABc(t^i,-c^i,°)^ l-e-/S"i (l-e-/^"0^ ' 

P , , i^3K,-^i,o) i^i'^(wi,-^i,o)e^"i 
/ci3A(o, -a;i,a;i) = ^ (e/3-i _ 1)2 ' 

and from Eqs. (|^3^ and lIX^ : 

? / N ^^2,3(^1,^2,^^3) J^i,2(wi,W2,tJ3) 

J^sc(a;i,a;2,a;3) = (^^^377:^)^^:^^ " (^^3.3 _ i)(i _ e'/^^O ' ^ ^ ^ 

r / ^ J?2,3(^l,^2,^3) ■*^1.2(W1,W2,W3) 

-'CBAlW3,a;2, ^ij — 



(l-e-'3"2)(i_e-/5'^3) (i_e-/3'^3)(e/3'^i _i)' 

that complete our main task and find spectral densities for multitime correlation functions from 
the known multitime temperature Green functions. 

4. Summary 

In conclusion, we have presented a general approach for derivation of the spectral relations for 
the multitime correlation functions. An analysis of the frequency dependences of their spectral den- 
sities is performed with special attention paid to the consideration of the non-ergodic (conserving) 
contributions. It is shown that such contributions can be treated in a rigorous way using multi- 
time temperature Green functions: representation of the Green functions in terms of the spectral 
densities and solution of the reverse problem — finding of the spectral densities from the known 
Green functions are given for the case of the three-time bosonic correlation functions. 

Generalization for the case of the higher-order multitime correlation functions and for the 
functions constructed also from the fermionic operators can be done in the same way, but it was 
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not be considered here. Here we can only note that for the higher-order functions, besides non- 
ergodic terms, which appear at zero value of one frequency and are connected with the conserving 
of one operator, the one connected with the conserving of products of operators, which appear at 
zero value of the sum of corresponding frequencies, always exist and they will be the only "non- 
ergodic" contributions for the pure fcrniionic correlation functions resulting from the presence of 
the matrix elements like {f\c^\l){l\c\f). Moreover, for some models the fermionic correlation and 
Green functions will contain only such "non-ergodic" contributions, e.g. local four-time two-electron 
Green's function for the Falicov-Kimball model [ i32. for which only contributions with zero sum 
of two frequencies exist. 

This publication is based on work supported by Award No. UKP2-2697~LV-06 of the U.S. 
Civilian Research & Development Foundation (CRDF). I am grateful to Professor I.V. Stasyuk for 
useful and stimulating discussions. 
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